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ABSTRACT 

White dwarf-white dwarf (WD-WD) and neutron star-neutron star (NS-NS) mergers may produce Type la 
supernovae and gamma-ray bursts (GRBs), respectively. A general problem is how to produce binaries with 
semi-major axes small enough to in significantly less than the Hubble time (tu), and thus accommodate the 
observation that these events closely follow episodes of star formation. I explore the possibility that such 
systems are not binaries at all, but actually coeval, or dynamical formed, triple systems. The tertiary induces 
Kozai oscillations in the inner binary, driving it to high eccentricity, and reducing its gravitational wave (GW) 
merger timescale. This effect significantly increases the allowed range of binary period P such that the merger 
time is fmerge < fH- In principle, Chandrasekhar-mass binaries with P ^ 300 days can merge in < fn if they 
contain a prograde solar-mass tertiary at high enough inclination. For retrograde tertiaries, the maximum P 
such that fnieige ^ IS yet larger. In contrast, P < 0.3 days is required in the absence of a tertiary. I discuss 
implications of these findings for the production of transients formed via compact object binary mergers. Based 
on the statistics of solar-type binaries, I argue that many such binaries should be in triple systems affected by the 
Kozai resonance. If true, expectations for the mHz GW signal from individual sources, the diffuse background, 
and the foreground for GW experiments like LISA are modified. This work motivates future studies of triples 
systems of A, B,and O stars, and new types of searches for WD-WD binaries in triple systems. 

Subject headings: stars: binaries: close, gravitational waves — celestial mechanics, stellar dynamics — 
stars: neutron, white dwarfs — supernovae: general 



1. EVITRODUCTION 

Stellar mass compact object mergers driven by gravitational 
wave (GW) radiation may power many types of observed as- 
trophysical transients and produce rare stellar exotica. No- 
tably, the merger of white dwarf-white dwarf (WD-WD) bi- 
naries has been suggested as a mechanism for producing Type 
la supernovae (Iben & Tutukov 1984; Webbink 1984), and 
neutron star-neutron star (NS-NS) or NS-black hole (BH) 
mergers are a leading model for the central engine of short- 
duration gamma-ray bursts (GRBs) (Ruffert & Janka 1999; 
Janka et al. 1999). Other types of compact object merg- 
ers such as those of BH-Helium star or BH/NS-WD binaries 
have been identified as possible mechanisms for long-duration 
GRBs (Fryer & Woosley 1998, Fryer et al. 1999ab), or pecu- 
liar supernovae (Metzger 2011). WD-WD mergers may also 
produce companionless millisecond pulsars (Saio & Nomoto 
1985) and R CrB stars (Webbink 1984, Nelemans et al. 2001). 
Finally, close compact object binaries brought near enough 
for mass transfer power a variety of astrophysical phenom- 
ena. One example is AM CVn stars (Warner 1995), which 
may produce faint ".la" supernovae (Bildsten et al. 2007). 

Estimates suggest that ^ 2% of the all the stars born with 
zero-age main sequence masses in the range ^ 2.5- 8 Mq 
become la supernovae (e.g., Horiuchi & Beacom 2010; but, 
see Maoz 2010). If sub-Chandrasekhar mass binaries can 
contribute to the la rate (e.g., van Kerkwijk et al. 2010), 
the IMF implies that a still smaller fraction of all stars pro- 
duce la's. Similarly, the massive star birthrate, which is ap- 
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proximately equal to the Type-II supernova rate, ^10^ Gpc~^ 
yr"^ (e.g., Horiuchi et al. 2009), is ~ 20-600 times larger 
than the NS-NS merger rate (Kalogera et al. 2001, 2004ab; 
O'Shaughnessy et al. 2008). 

These comparisons imply that although a large fraction of 
all stars are in binaries, only a small fraction have the req- 
uisite characteristics to produce compact object binaries that 
can merge via GW emission on a timescale short compared 
to the Hubble time. This issue is particularly important since 
there is evidence for a "prompt" subset of la supernovae that 
track the star formation rates of galaxies (Scannapieco & 
Bildsten 2005; Mannucci et al. 2006). "Prompt" may mean 

< 1 Gyr after a burst of star formation, or potentially even 

< 200 Myr (Aubourg et al. 2008; Brandt et al. 2010; Maoz 
& Badenes 2010; Maoz et al. 2010; Maoz et al. 2011). Sim- 
ilarly, short-duration GRBs occur in both quiescent and ac- 
tively star-forming galaxies (e.g., Berger 2009). In order to 
solve this problem, most studies focus on binary evolution 
channels that make more compact object binaries with small 
enough semi-major axis that a merger via GW radiation oc- 
curs "promptly" (e.g., Belczyhski & Kalogera 2001; Ruiter 
et al. 2009). Indeed, common envelope evolution of com- 
pact stellar binary systems can lead to compact object binaries 
whose GW merger time is significantly less than the Hubble 
time (Iben & Tutukov 1984, 1985, 1987). 

In this paper I forward an alternate hypothesis: the rare bi- 
nary systems that produce compact object mergers are not bi- 
naries at all, but instead hierarchical triple systems. In such 
systems, Kozai (1962) showed that the tertiary induces oscil- 
lations in the orbital eccentricity of the inner binary via a sec- 
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ular resonance. In the absence of general relativistic effects, 
the maximum eccentricity attained by the inner binary is 

eraax = ( l--COS^/j (1) 

where / is the incHnation of the outer orbit relative to the plane 
of the inner binary. Because the GW merger timescale is a 
very strong function of eccentricity (few « Pe- 
ters 1964, eq. |2l), the addition of a tertiary at high inclina- 
tion to a binary system can decrease few dramatically. Blaes 
et al. (2002) (hereafter BLS02) showed this explicitly for the 
case of triple systems of super-massive BHs, which would be 
formed by successive galaxy-galaxy mergers during structure 
formation. Miller & Hamilton (2002) independently applied 
this same idea to four-body interactions (binary-binary scat- 
tering) involving stellar mass BH binaries in globular clusters 
to accelerate the assembly of intermediate mass BHs in these 
systems. 

Here, I explore the possibility that WD-WD, NS-WD, and 
NS-NS mergers driven by GW radiation, of relevance particu- 
larly for la supernovae, GRBs, and other transients, are accel- 
erated by the presence of a tertiary at high inclination. I con- 
sider both coeval triple systems and triples formed by binary- 
single and binary-binary scatterings in dense stellar environ- 
ments. I show explicitly that the GW merger timescale for a 
subset of compact object binaries in triple systems, whether 
coeval or dynamically formed, is significantly decreased from 
the expectation for the binary alone. This effect dramatically 
increases the range of semi-major axes for which a merger 
will occur in a single Hubble time (fe)- Put another way, some 
systems which one might think have no hope of merging in tu, 
actually merge in <C Jh with a suitably placed tertiary. This 
will affect the population synthesis of la- and GRB-producing 
compact object binaries, and may significantly affect the over- 
all rate, even though there are fewer triple systems than lone 
binaries. 

These statements quickly raise the question of whether or 
not the parent population of triples that produce binary com- 
pact objects can accommodate the observed la or GRB rates. 
Although a complete discussion of this issue is beyond the 
scope of this paper, there are several reasons to believe that in 
the case of WD-WD mergers particularly, all such systems are 
(or were) triple. First, although only ^5 - 10% of all solar- 
type stars are thought to be in triple systems (Tokovinin et 
al. 2006; Raghavan et al. 2010), fully - 50% and - 100% 
of aU close solar-type binaries with orbital period P < 10 
and < 3 days, respectively, are in triple systems (Tokovinin 
et al. 2006; Pribulla & Rucinski 2006). Second, Fabrycky & 
Tremaine (2007) (hereafter FT07) have shown that the peak in 
the observed binary period distribution at ^ 3 days can be ac- 
counted for by the combined action of tidal friction between 
the two stars in the inner binary and Kozai oscillations in- 
duced by a hierarchical tertiary (see also Mazeh & Shaham 
1979; Wu & Murray 2003; Wu et al. 2007; Perets & Fab- 
rycky 2009). Thus, the closest binaries — those most likely 
to merge in less than a Hubble time — are precisely those that 
are most likely to be triples. Additionally, although little is 
known about the triple fraction of the A/B stars that give rise 
to WDs sufficiently massive to be plausible double-degenerate 
la progenitors and the massive O/B stars (> 10 Mq) that pro- 
duce NSs, there is evidence that the multiplicity of stars in- 
creases as a function of stellar mass (Lada 2006; Raghavan 
et al. 2010). Since the intermediate-mass and massive stars 



that produce WD and NS binaries become Cepheids during 
their post main-sequence evolution, the high occurrence of 
triples among such systems (Evans et al. 2005) further moti- 
vates consideration of their eventual compact objects.^ 

The high percentage of close solar-type binaries contain- 
ing tertiaries motivates this paper in part. The (albeit limited) 
statistics on such systems implies that the rate of la super- 
novae and GRBs compared to their overall parent population 
might be accommodated by the triple fraction (see 

There is an additional motivation for considering triple sys- 
tems as the progenitors for la supernovae. The delay-time 
distribution of la supernovae implies that many systems ex- 
plode on 1-10 Gyr timescales (e.g., Totani et al. 2008), and 
this fact implies a certain number density of progenitor sys- 
tems in the Galaxy. Simple estimates suggest that the nearest 
progenitor system is of order just ^ 50 pc away. So far, there 
have been no confirmed identifications of WD-WD binaries 
with combined mass greater than the Chandrasekhar mass dis- 
covered (see Mullally et al. 2009 for a recent compilation). 
Although there are a number of potential explanations for 
this fact — e.g., sub-Chandrasekhar mass mergers or single- 
degenerate systems dominate the la rate — an additional pos- 
sibility is that the reason no such systems have been found 
is that they are in triple systems. Since many such systems 
would be composed of a close WD-WD binary with a hier- 
archical main sequence tertiary, such systems would not have 
been found by searches that color-select WD binaries (Napi- 
wotzki et al. 2001; Badenes et al. 2009; Brown et al. 2010). 
Additionally, since the massive WD binaries that produce la 
supernovae may often have a lower-mass tertiary that evolves 
subsequently, triple WD systems may appear as a single WD 
(the young, low mass tertiary), and without a strong radial 
velocity signal from the distant, dimmer, but more massive 
WD-WD binary. Although a detailed discussion of these ef- 
fects on the selection of triple systems is beyond the scope of 
this work, I consider the lack of observed WD-WD progen- 
itors to be an additional motivation for exploring the strong 
triple hypothesis that all la progenitors are triple systems. 

The remainder of this paper is organized as follows. In ^ 
I discuss the types of systems that will be affected by Kozai 
oscillations. Although I outline some of the evolutionary pro- 
cesses that will affect triple-star evolution (e.g., common en- 
velope evolution and mass loss), for the purposes of this pa- 
per, I simply assume that binary compact objects are formed 
with a range of semi-major axes and masses, and with ter- 
tiaries with a variety of masses, semi-major axes, and inclina- 
tions relative to the inner binary. I then calculate the evolution 
of these representative systems using the methods described 
in ^and Appendix|A] The results are presented in ^ In ^ 
I discuss the results and the implications for transients pro- 
duced from these mergers, and the overall rate. I also discuss 
implications of these results for the gravitational wave fore- 
ground and background, as well as the signal from individual 
sources, for LISA (see also Gould 201 1). 

2. SCENARIOS 
2.1. WD-WD Binaries & Triples 

^ Typically, a Cepheid with a ~ 1 day period con'esponds to a ~ 3 Mq 
main sequence AO star, which produces a ~ 0.7 Mq WD (e.g., Kalirai et 
al. 2008), the minimum required for a Chandrasekhar mass, equal mass ratio 
binary. 
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Consider a system of two ^2-8 Mq (A5-B3) main se- 
quence stars (e.g., most commonly 2+2, 3+2 M©, etc., and 
perhaps preferentially "twins"; Pinsonneault & Stanek 2006), 
niQ and mi, in a close binary with semi-major axis ai, and a 
coeval hierarchical tertiary of mass 1112 and semi-major axis 
fl2. The minimum value of 02 is set by stability of the triple 
system: 02 <^ 3ai (Eggleton & Kiseleva 1995; Mardling & 
Aarseth 2001). 

If fli is less than a few AU, then the binary will undergo one 
or more mass transfer and common envelope (CE) evolution- 
ary phases, that decrease ai from its initial value by a factor of 
-10-100 (Iben & Tutukov 1984, 1985; Iben & Livio 1993; 
Ruiter et al. 2009). Although the physics of the CE phase(s) 
is uncertain, it is the primary mechanism for producing close 
WD-WD binaries that will merge in substantially less than tu 
(see eq.|2]l. 

In the triple system considered, mass loss from the inner bi- 
nary during the CE phase(s) will in general increase 02 from 
its initial value, typically by a factor of — 2 - 10, depending on 
the masses of system's constituents, and the final WD masses 
of the inner binary. Mass loss from m2 as it evolves to a WD 
will also increase 02- If 02 is less than a few AU either before 
or after the inner binary evolves, then there will be a com- 
plicated phase of triple CE evolution where the binary is en- 
gulfed by the tertiary, potentially decreasing both 02 and ai. 
A number of related evolutionary channels for triples were 
sketched by Iben & Tutukov (1999). 

In this paper, I show that for a range of parameters the 
timescale for the merger of WD-WD binaries via GWs is dra- 
matically decreased by the presence of a tertiary. As implied 
by equation ([T]i, and as discussed in detail in Sections|3]and|4l 
this mechanism requires ( 1 ) that the tertiary be at high relative 
inclination with respect to the inner binary and (2) that 02/01 
must be in the range of — 3 - 100. Both are expected to be 
modified by stellar evolution preceding WD formation. The 
latter is important since the combination of mass loss and a; 
evolution during the CE phase of a close binary will in general 
increase 02/01 from its initial value by a factor of 20- 1000. 
However, the phase of triple CE evolution that may occur as 
the tertiary evolves can in principle decrease 02 and oi, and in- 
crease or decrease 02/01. The evolution in inclination is also 
important since equation ([U and the fact that few oc (1 - e^)^/- 
imply that few will be a very strong function of /. For a ran- 
dom distribution of inclination angles between 0° and 90° 
(prograde tertiaries), the probability of having / > 70°, 80°, 
85°, and 89° is cos / ~ 0.3, 0.2, 0.09, and 0.02, respectively. 
Although the distribution of tertiary inclinations is unknown, 
the results of FT07 imply that if a uniform distribution in cos / 
is assumed when the stars first form, the resulting distribution, 
after many Kozai times (eq.|4|i, is quite flat (their Fig. 7), but 
with peaks near the critical Kozai angles (/ — 39° and 141°). 
Because tides detune the Kozai resonance, if the tertiary is 
at high /, the binary's semi-major axis and eccentricity will 
be affected by the combination of Kozai oscillations and tidal 
friction before WD formation (FT07). In particular, the inner 
binary will be driven to lower eccentricity than if Kozai acted 
alone (eq.[rii. If there is a CE evolutionary phase as the stars 
evolve, and mass transfer, the inner binary may well be nearly 
circular at the time of WD formation since the Kozai reso- 
nance will be terminated by such strong interaction between 
the inner components. However, subsequent to this compli- 
cated evolution, one expects Kozai to again operate after com- 



pact object formation. The system will then be driven to high 
eccentricity, and rapidly coalesce via GW radiation. 

Given the complications posed by triple star evolution, and 
in particular triple CE evolution, I simply assume that such 
systems lead to binary WDs with a range of a 1 , together with a 
tertiary WD or main sequence star with a range of inclinations 
and 02- The combination of triple star evolution with dynam- 
ics, mass transfer, and CE evolution, together with population 
synthesis is crucial from a theoretical perspective for evalu- 
ating the mechanism proposed here for rapid coalescence of 
WDs via gravity waves. However, these theoretical complica- 
tions aside, the results of Tokovinin et al. (2006) and Pribulla 
& Rucinski (2006) strongly motivate an observational search 
for WD-WD binaries with main sequence tertiaries. The dis- 
covery and characterization of such systems — in particular 
the quantities / and 02/01 — will be the ultimate arbiter in 
establishing the importance of the Kozai mechanism for rapid 
compact object mergers. 

Although I have made no attempt to model specific systems, 
recent work on the exotic planetary nebula SuWt 2 provides 
further motivation since it is thought to contain a close A- 
star binary with a (possibly WD) — 0.7 M0 tertiary (Exter et 
al. 2010). One can also imagine an important role for the 
Kozai mechanism in driving binaries composed of a WD and 
a main sequence star to contact, as in the single-degenerate 
scenario for la supernovae (Whelan & Iben 1973). 

Similar scenarios with different mass components may pro- 
duce AM CVn or R CrB stars, depending on the stability of 
mass transfer as the WDs interact (Webbink 1984, Nelemans 
et al. 2001), or millisecond pulsars (Saio & Nomoto 1985). 

In globular clusters, close WD-WD binaries may pick up 
a tertiary from the dense stellar field via either binary-single 
or binary-binary scattering. The latter should dominate the 
rate. This mechanism for the formation of triple systems 
that lead to Kozai-induced mass transfer has been studied 
by Ivanova (2008) and Ivanova et al. (2008,2010), and has 
been found to be important for producing the observed X-ray 
binary populations of globular clusters (see also Fregeau et 
al. 2004, 2009). As I show in ^ if the tertiary is captured 
into a high-inclination prograde orbit, or a retrograde orbit 
with / < 110- 120°, the GW merger timescale for the inner 
binary can be very short. Since such systems will have very 
high eccentricity, their GW signals will be peaked at high fre- 
quencies (®, and for individual systems LISA wiU see GW 
"pulses" at periastron (Gould 201 1). 

2.2. BH, NS-BH, NS, WD Scenorios 

For binaries with a BH or NS in triple systems, the su- 
pernova explosions that accompany NS (and possibly BH) 
formation may in some cases unbind the inner binary or 
the tertiary, depending on the binary mass ratio, the tertiary 
mass, the eccentricities of the inner and outer orbits, and 
the magnitude and direction of the "kick" velocity given to 
the NS at birth. For NS-NS-producing binaries of relevance 
for short-duration GRBs, one imagines an inner binary with 
components of -9-12Mq (e.g, 9+10, 9+9, II+IOMq, 
etc.), forming via standard scenarios (Bhattacharya & van 
den Heuvel 1991), and with a relatively distant tertiary with 
m2 — 5 - 8 M0, which may be either a main-sequence B star, 
or perhaps a BH formed by a more massive star Depend- 
ing on the mass of the system expelled before and during the 
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supernovae that produce the NSs and the magnitude and di- 
rection of their kicks, the tertiary may remain bound, but with 
small fl2 preferred (see, e.g.. Fig. 10 of Kalogera 1996 for the 
binary case; Hills 1983). 

After the supernovae, the resulting NS-NS binary will be 
subject to Kozai oscillations by the tertiary, depending on its 
semi-major axis and inclination. If the tertiary is a main- 
sequence star, it will become a massive WD with ^ 1 - 
1.4M0. As in the WD- WD case, because the timescale for 
mass loss during this transformation is long compared to the 
orbital times in the system, it should remain bound, albeit with 
larger a2/a\. In addition, the system may also undergo triple 
CE evolution as the tertiary evolves. After this is complete, if 
the system persists, Kozai cycles will then resume. 

Many similar scenarios can be considered for NS-WD bina- 
ries or BH-NSAVD binaries in coeval triple systems, although 
the survival of the triple through the formation of the compact 
object binary must be evaluated case-by-case. A discussion 
of one such system in the context of Kozai oscillations in a 
hierarchical triple system is given in Champion et al. (2008) 
and Freire et al. (20 11) for the system PSR J1903H-0327. 

As discussed in g2j]in the WD-WD case, NS-NS/WD/BH 
binaries that undergo Kozai oscillations can also be formed 
dynamically in dense stellar environments via binary-single 
and binary-binary scattering. A significant fraction of such in- 
teractions would be expected to produce a stable triple system 
(see MH02; Wen 2003; Ivanova 2008; Ivanova et al. 2008, 
2010). 

3. METHOD 

3.1. Timescales 

I consider an inner compact object binary with semi-major 
axis a\, eccentricity ei, argument of periastron^i, and masses 
m^) and m\ . The hierarchical tertiary has mass m2, semi-major 
axis a2, argument of periastron g2, and mutual inclination with 
respect to the inner binary of /. Subscripts of "0" are used to 
indicate initial values (e.g., ei o). 

There is a strong hierarchy of timescales in the problem of 
triple systems consisting of a compact object binary. First, in 
the limit of high eccentricity, the GW merger timescale of the 
inner binary is (Peters 1964) 
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where M = mQ + m\. The low-eccentricity version of equation 
(H is few X (425/768)(l-e2)-7/2 (Peters 1964). Through- 
out this paper I focus on stellar-mass binaries whose initial 
semi-major axis is large enough that the nominal value of 
few, in the absence of a tertiary, is greater than the Hubble 
time, t\i ~ 14Gyr This implies values for oi o larger than 
^ 0.017 AU, depending on the masses of the binary compo- 
nents considered. 

As discussed by BLS02 and MH02 in the context of super- 
massive BH mergers and the formation of intermediate mass 
BHs, respectively, GR periastron precession "de-tunes" the 
secular Kozai resonance. This effect in general decreases the 
maximum eccentricity attainable at fixed tertiary inclination 



(see Appendix of BLS02 and the discussion in FT07). The 
timescale (period) for GR precession (e.g., eq. 23 of FT07) is 
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and the Kozai timescale is (Innanen et al. 1997; Holman et 
al. 1997) 
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where b2 = a2{l — e'^^l^. As discussed in BLS02, Kozai oscil- 
lations only operate if fK < fcRp- The strong dependence of fK 
on fl2/fli implies that there is a maximum fl2, beyond which 
Kozai oscillations are ineffective. Additionally, since 

fCRp/fK OC fli, 

as the binary evolves to smaller fli because of GW radiation, 
fK eventually becomes larger than fcRp- This can cause the 
binary to circularize before coalescence, and it is the compe- 
tition between fcRp and fK that determines much of the time 
evolution of the system as it evolves towards merger 

Momentarily ignoring the complications of GR precession, 
one can combine equation ([T]i with equation ^ to get a rough 
order-of-magnitude sense of the importance of Kozai oscil- 
lations in triple systems for the rapid merger of compact ob- 
jects. Taking the maximum eccentricity the system reaches to 
be gmax, the roughest estimate of the merger time is simply 
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where the factor (1 -Smax)"'''^ corrects for the small relative 
amount of time the system spends at high eccentricity (e.g., 
MH02). Substituting, one finds that 
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which shows the very strong expected dependence on the in- 
clination angle. The estimate of equation ^ is only valid 
for angles in the critical Kozai range between 39° < / < 141°, 
and fails to account for important corrections to e,nax from GR 
precession (BLS02, MH02, Wen 2003, FT07). In particular, 
it has no dependence on 02/01, and its dependence on fli is 
only approximate. In some regions of parameter space equa- 
tion (|6]l thus grossly underestimates fmerge (see ®. A much 
more accurate, but less simply stated, estimate of fmerge can 
be made using the method of Wen (2003) that is accurate to a 
factor of a few over many decades in fmerge (see AppendixlAb. 

3.2. Equotions, Assumptions, & The Merger Time 

I solve the octopole-order equations for the secular evolu- 
tion of the orbital elements of the system, as given in BLS02 
(their eqs. 11-17), which are based on the expressions derived 
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Fig. 1 . — Time evolution of a triple system with fli.o = 0.05 AU and aa.o/'^i.o = 20, with initial conditions (q = 85°, ei.o = ^2.0 = 0.1, gi o = 0°, g2.o = 90°, 
niQ = 0.8 M0, m[ = 0.6 Mq, and 013 = 1.0 Mq (i.e., a WD-WD binary with a main sequence companion). The nominal GW merger timescale without the tertiary 
is ~ 3 X 10'^ yr, but with the tertiary, Wrge — 2.5 X 10* yr. Left panel: cos(i) (dotted), cos(gi) (dashed), ej (solid). Right panel: e[ (black), oi/oiq (gray). 




lO'o - 



109 ^ 



10= 



10' 



10" 



10^ 



10* 



10= 





1 1 1 1 1 1 








lGv,{a,,o=0 01 AU, e,„=0.1) n 








































\ ^\ 40^^ 






\ 


/ / / B 
// / 
/ / / V<> 






















\ ^ \ 


/ ' s!i.o = 0.1 e3„ = 0.1 ! 






/ / gi.o = 0" 83,0 = 90" : 




\ \ 
\ \ 


/ / 

/ ' m„=0,aMg E 




\ \ 


/' m, = 0,6Mo - 






/ m,= l,OM^ ^ 

/ : 




\ .1/ 


i - 

1 - 



-0.4 



-0.2 



^2,0/^1 





cos(i„) 



0.2 



0.4 



Fig. 2. — Left Panel: Merger time fmcrge versus initial 02.0/01.0, for oi.o = 0.01 (light solid), 0.1 (dashed), and 0.5 AU (heavy solid), for several values of (q. 
For sufficiently large 02,0/01,0, the presence of the tertiary does not speed up the merger and fmergc — ^ fGw(0| Oi^l.o = 0.1). Right Panel: fmcrge versus cos(io) 
for oi = 0.01 AU, and 02^0/01,0 = 20, 30, 40, 50, 60, 70. Note that for retrograde tertiaries Wrge becomes very short, and very high eccentricity WD-WD 
interactions are possible. The scaling of /merge oc cos* (' from equation j6) is shown for comparison. 



in Ford et al. (2000, 2004) (Krymolowski & Mazeh 1999; 
Marchal 1990). These equations include both GR periastron 
precession and GW radiation (MH02; Wen 2003), but ne- 
glect tidal forces and treat the masses as point particles. This 
amounts in part to neglecting terms in the equation for the 
time evolution of the longitude of periastron of the inner bi- 
nary associated with both tidal and rotational bulges. Both can 
suppress Kozai oscillations in a manner similar to GR preces- 
sion (eq. |3]i. In particular, the timescale associated with the 



apsidal motion induced by a tidal bulge is (e.g., FT07) 
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where R is the radius of the compact objects (assumed equal), 
mo = mi has been assumed, and k is the classical apsidal mo- 
tion constant (for the numerical estimate A: = 0.1). For suffi- 
ciently high eccentricity (e.g., ei > 0.9986 for the parameters 
of eq. IT), ftide becomes less than fcRp, and one expects tides 
to become important to the evolution. In some cases the ec- 
centricity can become large enough that periapsis approaches 
the WD radius. In such cases, one expects strong tidal heat- 
ing and circularization. Although these effects are of interest 
in their own right, and they might be especially important for 
the resulting transients from such mergers, they are not cap- 
tured in the current study, and will be the subject of a future 
work. For the present purposes, it is sufficient to note that 
GW merger timescale on a scale comparable to the compact 
object radius (whether NS or WD), is ultrashort compared to 
the merger time of circular binaries without a tertiary at the 
semi-major axes of interest (> fe). Thus, if the binary does 
circularize on such scales, the merger timescale fmei-ge as esti- 
mated by numerical solution of the time-dependent equations, 
is not dramatically affected. 

While the merger time fmerge reported for binaries calculated 
in mis the time required for the semi-major axis to reach the 
radius of the compact object Rqo (I take Rwd = 5000 km and 
Rns = 10 km), in some extreme cases the periapsis of the orbit 
reaches ^ Rco- In such cases, I explicitly note that such solu- 
tions will be strongly affected by tides. The combined action 
of tidal friction and Kozai oscillations on the period distri- 
bution of solar-type binaries in triple systems has been con- 
sidered in detail by FT07 (see also Mazeh & Shaham 1979; 
Kiselevaet al. 1998; Eggleton & Kiseleva-Eggleton2001; Wu 
& Murray 2003; Wu et al. 2007; Perets & Fabrycky 2009). A 
similar calculation of white dwarf and neutron star binaries in 
triple systems is saved for a future paper. 

Finally, all of the initial configurations described below are 
chosen to be stable when compared to the empirical 3-body 
stability criterion of Mardling & Aarseth (2001) (see also 
Eggleton & Kiseleva 1995). For a discussion in the context 
of Kozai oscillations, see BLS02. 

3.3. Numerics 

Standard methods (e.g., Bulrisch-Stoer; Press et al. 1992) 
are employed to solve the system of equations, and I have 
done several checks to ensure the fidelity of the results pre- 
sented. First, I have varied the numerical tolerance of the al- 
gorithm systematically and found that the results presented 
here are converged. As an additional check, I have verified 
that as fl2 becomes large and fK ^ fcRp, the solution for f merge 
of the inner binary approaches the result of Peters (1964) for 
an isolated binary. Finally, I have spot checked my calcula- 
tions directly against the code of BLS02 and the numerical 
results of Wen (2003), and find excellent agreement. 

Some of the results presented make use of the simple ap- 
proximate method of Wen (2003), described in Appendix lAl 
Comparisons between this approximation and the actual solu- 
tion of the time-dependent problem are provided there. 

4. RESULTS 

The parameter space of possible masses and orbits is very 
large. To restrict the total model space, I take the initial val- 
ues of the orbital eccentricities and arguments of periastron to 
be ei_o = 62,0 = 0.1, gi,o = 0°, and §2,0 = 90° throughout this 



paper. One expects the assumption of low initial eccentricity 
to be reasonable except in two cases of particular interest: (1) 
in dynamically formed triple systems the eccentricity of the 
tertiary may be large 62 ^ 0.9 (Ivanova 2008), and (2) in NS- 
NS binaries the NS kicks at buth may cause ei to be large. 
Higher e\ and 62 generically lead to faster mergers, and thus 
the assumption of low initial eccentricities is conservative. 

Figure [T] shows the time evolution of a selection of orbital 
elements for the fiducial WD-WD case with mo = 0.8 M0, 
mi = 0.6 M0, and m2 = 1 .0 M©. The masses of mo and mi are 
chosen to sum to 1 .4 M0 for illustrative purposes, and mo T^^mi 
so that the octopole-order terms in the dynamical equations 
operate. Results for mo = mi = 0.7 M0 are not significantly 
different. The initial semi-major axis (fli.o = 0.05 AU) is cho- 
sen so that the nominal GW merger timescale in the absence 
of the tertiary is ^ 200fH- The initial mutual inclination of 
the system is / = 85° (cos / ~ 0.09). The left panel shows the 
early time evolution of / (dotted), g\ (dashed), and e\ (solid). 
The periodic changes in e\ correspond to Jk (eq.lUi. The right 
panel shows the late-time evolution of ex and the semi-major 
axis, scaled by its initial value (ai/ai,o). The system merges 
in ~ 2.6 X 10^ yr, approximately lO'* times faster than without 
the tertiary. At no time before the very end of the calculation 
does the periapsis of the orbit reach 2 x Rsnd- 

The results of many such calculations are presented in Fig- 
ure |2l The left panel shows fmeige as a function of the ini- 
tial value of fl2.o/fli.o, for oi o = 0.01 AU (thin solid), 0.1 AU 
(dashed), and 0.5 AU (heavy solid), for several different val- 
ues of /o as labeled. The nominal (no tertiary) binary merger 
timescales via GW radiation are ~ 5 x 10^ yr (^ 0.3 fn), ~ 
5 X 10'^ yr (3 x IO^Jh), and - 3 x 10'*'yr (2 x 10^ fe), for 
fli.o = 0.01 AU, 0.1 AU, and 0.5 AU, respectively. 

In the models with oi.o = 0.01 AU, as fl2.o/fli.o becomes 
greater than ~ 50-60, there is essentially no decrease in fmerge 
with respect to the case without the tertiary, and all models ap- 
proach fnierge ~ 5 X lO'^ yr. For fl2,o/ai.o ^ 20, there is a mini- 
mum in fmeige- The retrograde cases shown with / = 95° have 
very short fmerge for small fl2,o/ai,(), and typically merge in a 
single Kozai timescale Jk. 

Setting the WD radius to be /?wd = 5000 km, for the cases 
with a2,()/fli,o = 10 and 20 and / = 89°, the periapsis of the in- 
ner binary orbit becomes less than the 7?wd in the first Kozai 
oscillation, at time fK- However, for 7?wd = 1000km it does 
not, and the evolution is qualitatively similar to that presented 
in the right panel of Figure [T] Clearly, for these cases a more 
complete model with tidal dissipation and circularization is 
required to capture the dynamics and to make an accurate cal- 
culation of fmeige- For the purposes of constructing this figure, 
I have assumed that /?wd is small enough that a "collision" 
(periapsis < /?wd) does not occur, and thus the results pre- 
sented may be an upper limit to fmerge - 

The right panel of Figure |2] shows fmerge as a function of 
cos/ for fli.o = 0.01 AU and 02.0/01,0 = 20, 30, 40, 50, 60, and 
70. The estimate of equation (|6l), which fails to capture the 
very strong dependence on 02/01 is shown as the dashed line. 
Again, for 02,0/01,0 < 20 and 89 < /o < 96, strong tidal in- 
teractions in a single fx are expected; for these models, as e\ 
reaches its first maximum, the periapsis is less than the fidu- 
cial WD radius of 5000 km. However, as in the left panel, the 
results shown assume /?wd small enough that a "collision" 
never occurs. All other models have the same qualitative be- 
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Fig. 3. — Range of allowed initial inner binary semi-major axis ai versus initial mutual inclination io, such that the inner binary merges in ?h, 'h/10, Jh/IOO, 
and ;h/1000 (darkest to lightest), for outer tertiary semi-major axis «2.o/'3i.o = 10 (top left), 40 (top right), 60 (bottom left), and 80 (bottom light), computed by 
estimating fmcigc using the algorithm described in Appendix |A]for mo = mi = 0.7, 1112 = 1.0, ei o = £2,0 = 0.1, gi o = 0°, and ^2.0 = 90°. Note that for these inner 
binary masses, only those with a 10 < 0.015 AU (off the bottom of the range shown) would be expected to merge in ?h in the absence of a tertiary companion. 
Note further that for tertiaries randomly distributed in (', cos(0 is proportional to the probability of having such a system. 

0.7, and mi = 1 .0. As a function of the initial value of the in- 
ner binary semi-major axis, oi.o, and as a function of cos/q, 1 
survey this parameter space for regions where fmerge is < fu 
(black), < O.lfH (dark grey), < O.OlfH (hght grey), and < 
O.OOlfH (interior white) for initial values of fl2,o/fli,o = 20 (up- 
per left), 40 (upper right), 60 (lower left), and 80 (lower right). 
Only regimes where the nominal GW merger timescale with- 
out the tertiary is larger than tu are explored (ai o > 0.015 AU; 
see eq. |2]l. The right vertical axis gives the associated no- 
tertiary GW merger time; considered values for ai o run up to 
lOAU, equivalent to few ^ 10^^ yr without the tertiary. The 
basic trends presented in the dynamical calculations shown in 
Figure |2] are reproduced. Exceedingly rapid mergers can be 
induced for a narrow range of retrograde values of / near 95°, 



havior as shown in Figure [T] 

In order to make a broad exploration of parameter space 
for many models, instead of calculating the detailed time evo- 
lution of each system, as in Figures [T] and |2] 1 use the ap- 
proximate method described in Wen (2003), which 1 detail in 
AppendixlAl(see discussion after eq.|6]l. As shown in the Ap- 
pendix, the method generally underpredicts fmerge by a factor 
of ^ 1.5-2 in most regions of parameter space, but it can 
overestimate fmerge by a factor of ^ 10 for some retrograde 
cases. In general, the method is simple, robust, and accurate 
over many decades in f merge • 

Using this approximate method. Figure [3] shows results 
analogous to those presented in Figure |2l but for niQ = nii = 
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Fig. 4. — Total mass of compact object binaries M = mo versus orbital period (lower axis) and semi-major axis (upper axis, assuming M = 1.4Mq). The 
heavy solid lines show the critical value of the inner binary period (P™"") such that fmcigc = 'h, including a tertiary with iq = 70, 80, 85, 89, and 90 degrees, 
assuming mo = mi, 1112 = I.OMq, and the same parameters used in Fig. [3] Each line is labeled with the values of the critical tertiary period P2 (yrs) for which 
pmax occurs (small numbers). Dotted lines are of constant no-tertiary GW merger time of ?gw = 1 (heaviest), 10^^, 10*, and lO' fH, assuming mo = mi, and e[ = 0. 
Thus, for example, triple systems consisting of binaries with M = 1.4M0 (along the horizontal dashed line) and Pi < 500 days can in principle merge in ?h- For 
larger 1112, the allowed range of Pi increases. For some retrograde orbits (io ~ 95°), it increases dramatically (see Fig. [3). Data on WD-hWD (filled squares), 
NS-l-WD (filled-Hopen squares), and NS-l-NS (fiUed-Hopen circles) binaries are shown (Mullally et al. 2009; Kulkarni & van Kerkwijk 2010; Nelemans et al. 2005; 
Kilic et al. 2010ab, Stairs 2004). The Chandrasekhar mass (dashed line) and the maximum NS mass (grey shaded) are indicated for reference. 



and for many prograde tertiary orbits. As 02,0/01,0 increases, 
the allowed area in the ai.o-cos/o plane decreases, and for 
02.0 /fl 1.0 Jc^ 100, Kozai oscillations do not accelerate mergers; 
this follows from the strong dependence of on 02/01 and 
the effects of GR precession. Nevertheless, it is clear from 
Figure|3]that for a significant region of parameter space, fmerge 
can be less than 10, 1, 0. 1, or even 0.01 Gyr, even for systems 
that have nominal (no tertiary) merger times of fe- 

Using these estimates of fmerge, it is possible to calculate the 
maximum possible value of the inner binary semi-major axis 
o™" such that merger occurs in a single Hubble time, given a 
tertiary of mass 1112 =1.0 Mq at any /. Figure[3]implies that for 
retrograde orbits o"™ is extremely large near iq ~' 95°. If we 
restrict our attention to prograde orbits, there is a unique value 
of o™" for each initial io such that fmerge < fn- As implied by 
the left panel of Figure |2] since fmerge exhibits a minimum as 



a function of 02.0/01.0, o™", will occur not at the smallest 
02.0/01,0, but instead at an intermediate value determined by 
the shape of the minimum in f merge (02.0/0 1.0) ■ 

In Figure m the solid lines show the critical maximum inner 
binary period P™" corresponding to this o™^" as a function of 
M = mo + nil, for tertiary mass of 1112 = l.OM© and inclina- 
tion of 70°, 80°, 85°, 89°, and 90°. Period is shown instead 
of semi-major axis to make contact with observational plots 
of compact binaries prevalent in the literature. The top axis 
gives o 1 for M = 1 .4 Mq for reference. The dotted lines are of 
constant GW merger timescale assuming ei = and mo = mi 
for multiples of fe- In the absence of a tertiary, forM = 1 .4 M© 
(dashed line) a binary period of < 0.4 days is required for 
merger in tu- The solid lines show that the allowed range 
in observed inner binary period is dramatically increased for 
binaries with a hierarchical tertiary. The small numbers along 
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the solid lines denote the value of the tertiary period P2 (using 
M+m2) at each value of f™" (a™") for each /. Again taking 
a Chandrasekhar mass inner binary with M = 1 .4 (follow- 



ing the dashed line), for / = 70°, I find that a" 



: 0.037 AU, 



pmax 
^1 



2.2 days, and aj/ai ~ 37, implying that the 1 Mq ter- 
tiary is required to be at a2 — 1.37AU, with an orbital pe- 
riod of P2 ~ 380days (using M + mi = 2.4M0). As another 



example, for / = 85°, a 



max 
1 



0.175 AU, A" 



; 22.7 days, 

a2/fli ~46, fl2 ~ 8AU, andf2 — 5400days ~ 15 years. 

Figure |4] implies that for any given observed WD-WD bi- 
nary with M ~ 1 .4 Mq, in order to exclude the possibility that 
the system will merge within one Hubble time, one must ex- 
clude the possibility of a tertiary with some mass 1112, at some 
inclination /. Given the very large values of P2 for / > 80°, 
this may be observationally challenging. Almost all of the 
WD-WD binaries could in principle merge in less than fn if 
one posits the existence of a ^ 1 Mq tertiary at > 80° and with 
02/01 — 40 (P2 ~ 4 yr). Note that the values for f™" shown in 
Figure m shift to yet longer periods if 1112 increases. For retro- 
grade orbits, particularly near ^95° (see Fig. [3]), they increase 
still further Again, some caution is warranted in interpreting 
the solid lines in Figure|4]for / > 89° since for WD-WD bina- 
ries tidal interactions will almost certainly affect the estimate 
of fmerge- Nevertheless, it is difficult to see how tides would 
dramatically increase fmerge when the relevant comparison is 
torn. 

4.1. NS-NS Binaries 

Similar plots to Figure [3] may be generated for the case 
of NS-NS mergers. Assuming mo = m\ = 1.4Mq, the left 
panel of Figure |5] shows the region in the oi q-cos/q plane 
for which a merger occurs in ^ 5 x lO'^yr ~ fH/300, which 
approximates the main-sequence lifetime of the tertiary, as- 
sumed here to be m2 = 6 Mq . The right panel shows the same 
quantities, but for m2 = 1 Mq , as might be appropriate after the 
tertiary becomes a WD. As discussed in ^2.21 in some cases 
the tertiary might also be a massive BH, in which case the al- 



lowed region for which fmerge is less than fn grows. Similar 
to the WD-WD case, Figureg] shows that for a 1 Mq tertiary 
with / > 85° and P2 < 10 yrs, all the observed NS-NS systems 
merge in fn given an appropriately placed tertiary. 

The parameter space for which the tertiary remains bound 
to the system through both of the supernovae that produce the 
NS-NS binary may not be large. In particular, the tertiary can 
neither be too close, because of dynamical stability, nor too 
far way because of the binding energy of the orbit relative to 
the kick imparted to the center-of-mass during each of the in- 
ner binary component's supernovae. However, inspection of 
the plots by Kalogera (1996) for the binary case indicate that 
for relatively small tertiary semi-major axis (02 < few AU), 
an intermediate-mass tertiary can stay bound, and perhaps 
acquire high eccentricity during the supernovae for kick ve- 
locities of ^ 100km s"' in some cases. A high value of 62 
could lead to more rapid mergers than implied by Figure |5] 
As discussed in Sections [2.11 and l272l 02 ^ AU likely guaran- 
tees strong (common envelope) interaction between all com- 
ponents as the primary and then the secondary evolve off the 
main sequence. 

5. DISCUSSION 
5.1. Rates 

The fact that essentially all close solar-type binaries are in 
triple systems argues that their subsequent WD-WD binaries 
will also be in triple systems. Depending on the semi-major 
axis and inclination distribution of the tertiaries at the time 
of formation of the compact objects one expects Kozai os- 
cillations to speed up the process of coalescence of the inner 
binary significantly, as shown in Figures [T][3] The work of 
FT07, together with Figure[3] serve as a guide to an estimation 
of the rate, but the estimate is complicated by the expectation 
that the distribution of inclinations may be biased towards co- 
planar orbits for coeval systems. 

Nevertheless, to make a simple rate estimate, assume that 
the distribution of binary and tertiary semi-major axes is flat 
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(equal numbers in log semi-major axis), and that the probabil- 
ity of having a system at inclination / is dp /di = sin /. Impor- 
tantly, Figure[3]shows that the range of oi strongly affected by 
Kozai oscillations is comparable to the relevant range of close 
binary semi-major axes. In addition, the Kozai mechanism 
operates over ~one decade in 02/01 ~ 3 - 100, which corre- 
sponds to ^ 2 decades in 02- Thus, one then expects many of 
the triples to be strongly affected. Finally, although the area in 
the fli - cos / plane is smaller for more rapid mergers, the sys- 
tems at large ; dominate the rate. For example, comparing the 
prograde regions of the upper left panel of Figure [3] one sees 
that the total area of the shaded regions (with ai measured 
in log units) measures the fraction of all systems with that 
merger time. Although the area of the black region is roughly 
10 times that of the white, the latter merge 1000 times faster 
Thus, the average rate is dominated by the highest inclina- 
tion systems, a consequence in part of the very steep scaling 
of Wige with cos/ (e.g., eq. |6l Fig. |2]i. This can be shown 
explicitly by using the very crude estimate of equation (|6]l. 
Momentarily ignoring the dependence on 02, one may write 
the total merger rate as 
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where is the total number of systems. As in equation (|6]l, 
this expression neglects the very important 1112 dependence of 
the Kozai mechanism and the 02/01 dependence, and thus 
equation dHJ should be considered schematic. For practical 
purposes the dependence on oa/fli can be considered nearly 
a step function since for all relevant parameters Kozai is inef- 
fective for 02/01 > 100, as shown in Figures|2]and[3] As dis- 
cussed in Section |2] the latter is particularly important since 
the inner binary will likely undergo two common envelope 
events and significant mass loss, thus increasing 02/01 from 
its initial value. Nevertheless, if triple systems are prevalent, 
one expects them to contribute both to the very prompt rate 



(< lO'^yr; 



and to the delayed rate (>Gyr), since for 



• < 39 2° (and/or 02/01 > 100 and/or 1112 < 1 Mq), Kozai 
does not affect the merger time of the inner binary. Thus, 
for low-/ or large 02 /o[ systems, one should obtain a delay- 
time distribution identical to the binary-only case. A more 
complete calculation of the delay-time distribution of triple 
systems is saved for a future work. 

5.2. How "prompt" is prompt? 

A prompt component to the supernova la rate has been 
claimed by a number of authors (e.g., Scannapieco & Bild- 
sten 2005; Mannucci et al. 2006; Aubourg et al. 2008; Brandt 
et al. 2010; Maoz & Badenes 2010; Maoz et al. 2010; Maoz 
et al. 2011). Various estimates suggest that ^ 50% of la's 
are "prompt," with a characteristic short timescale of ~' 0.1 - 
1 Gyr, and potentially with a two component power-law or bi- 
modal delay-time distribution (see Mannucci et al. 2006 and 
the recent review by Maoz 2010). 

In the picture presented here, a fraction of intermediate- 
mass stars are born in triple systems. The fraction with high 
inclination tertiaries can merge extremely rapidly as soon as 
both stars in the inner binary are WDs. Because the calcu- 
lated merger times are in many cases much less than even 
the post-main-sequence timescales of these stars, one expects 
the fastest WD-WD mergers to be limited only by the main- 
sequence lifetime of their progenitors. Thus, ^ 8 + 8 M0 bina- 
ries in triple systems merge first, immediately after WD birth. 



and thus the minimum delay time is ^ 3 x 10^ yr. In this case, 
one expects a delay-time distribution of ^ f"'/^ (Pritchet et 
al. 2008). These prompt supernovae would preferentially be 
super-Chandrasekhar mass binaries (see Pakmoret al. 2010). 

Many of these statements are equally applicable to NS- 
NS mergers. For example. Figure |5] shows that the merger 
timescale for NS-NS binaries in triple systems can be very 
short. Again, there is a region of parameter space where the 
limiting factor is the time required to produce the NS-NS bi- 
nary, again implying that the fastest NS-NS mergers can come 
just ^ 10^ yr after the last star formation episode. This may 
explain the fact that many short-duration GRBs are seen in 
star-forming galaxies (e.g., Berger 2009). By extension, one 
expects similar arguments to hold for BH-WD/NS or NS-WD 
systems that in some cases might give rise to long-duration 
GRBs (Fryer et al. 1999). 

These considerations are amplified if the multiplicity of 
stars increases with zero-age main sequence mass (Lada 2006; 
Raghavan et al. 2010). 

5.3. Progenitors 

Even if they constitute much of overall merger rate, the pro- 
genitor systems for the shortest-lived compact object mergers 
will be very rare in a given galaxy simply because the time 
a given system spends as a progenitor is very short. Never- 
theless, as in single-degenerate models of la supernovae, for 
WD-WD systems with Gyr merger time (with the tertiary!), 
one could look for close binaries with evidence of a tertiary 
companion. This would work if the tertiary is a WD or NS, 
but the problem is that in many cases one expects the com- 
panion to be a main-sequence star of ~ 1 M© (® . The latter 
would necessitate a new search strategy, since most searches 
for close WD-WD binaries are color-selected (Napiwotzki et 
al. 2001; Badenes et al. 2009; Brown et al. 2010). One exam- 
ple would be to target intermediate mass main sequence stars 
for radial velocity measurements and search for the signal of a 
more massive, but unseen companion, which might be an old, 
high-mass WD-WD (or NS/BH) binary. 

The same is true of the coeval NS-NS systems described in 
^ In that case, one might look for NS-NS binaries in pulsar 
searches, with timing characteristics that suggest the presence 
of a tertiary component. However, in many cases the NS com- 
ponents may not be observed as pulsars, and the tertiary may 
be a main sequence ^6-8 M© star. In both the WD-WD 
and NS-NS cases, the only visible progenitor at the sight of 
the subsequent explosion (the compact object merger) in pre- 
explosion imaging may be a bright main sequence star (sim- 
ilar in spirit and conclusion to the recent work by Kochanek 
2009). 

5.4. Tronsients 

Recently, several authors have discussed the possibility that 
WD-WD collisions, at small impact parameter, may be a way 
of producing la-like supernovae (e.g., Rosswog et al. 2009; 
Raskin et al. 2010). For some of the calculations presented in 
Figures|2]and|3]the merger occurs in a single Kozai timescale, 
at very high eccentricity. I have neglected tidal forces on the 
compact objects throughout this work (® , but these findings 
suggest that in some cases something akin to a "collision," or 
at least a very strong tidal interaction, may be induced by the 
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Kozai mechanism. This is particularly promising for retro- 
grade orbits with / ^ 95° (Fig.O. Such tertiaries may be cap- 
tured in binary-binary collisions in dense stellar environments 
(Ivanova et al. 2008, 2010), or perhaps in some cases may be 
coeval (see Fig. 7 of FT07). Simple numerical experiments 
with the evolutionary equations described in ^ but including 
apsidal motion of the inner binary as a result of tides (eq. U}, 
still allow for a region of parameter space where the merger 
occurs at very high eccentricity. 

5.5. Interactions & The Tertiary 

WD-WD mergers: If WD-WD mergers driven by a hier- 
archical tertiary produce Type la supernovae, the explosion 
should interact with and overtake the tertiary. Thus, just as 
in searches for the remaining star in the single-degenerate 
scenario (e.g., Kerzendorf et al. 2009), in the triple scenario 
proposed here, the tertiary in historical Galactic supernovae 
should be currently inside the la supernova remnant. In addi- 
tion, depending on the structure of the tertiary, some amount 
of mass might be expected to be lost by the Shockwave in- 
teraction (Marietta et al. 2000). However, in contrast to the 
single-degenerate picture (Whelan & Iben 1973), one expects 
the tertiary to be in many cases more than an AU distant from 
the explosion and would only very rarely be a giant. Us- 
ing the simple analytic model of mass stripping and ablation 
in supernova explosions with binary companions developed 
by Wheeler et al. (1975), 1 find that even in the most com- 
pact cases considered here, with 02 ^ 0.1 AU, the total mass 
stripped and ablated from a main sequence solar mass tertiary 
is just ^ 0.01 M0. Forfl2 1 AU, the expected mass lost from 
the tertiary will be minimal.^ Similarly, based on the work of 
Kasen (2010), one expects a soft X-ray flash as the la shock- 
wave overtakes the tertiary on a timescale of ^ 10"* s ((22/ AU) 
for a Shockwave of ^ 10'* km s~' and a very small change to 
the early-time optical and UV lightcurve. 

NS-NS mergers: If NS-NS mergers in triple systems pro- 
duce short-duration gamma-ray bursts (SGRBs), then there 
is a chance for interaction of the relativistic blastwave with 
the tertiary companion because the systems most likely to 
be affected by Kozai oscillations are at high mutual incli- 
nation relative to the inner binary, and because the inferred 
opening angle for SGRBs is relatively large (e.g., 10-30°; 
Nakar 2007). One then envisions an interaction similar to 
that calculated in MacFadyen et al. (2005), but with either 
a main-sequence, WD, or BH companion to the SGRB ex- 
plosion (see ^2.21 Fig. |5]l. The timescale for interaction is 
^ aijc ^ 500(fl2/AU)s. Moreover, because we are presum- 
ably looking roughly down the jet axis, the tertiary may in 
some cases be roughly along the line of sight. In the case of 
a main-sequence companion, some material may be stripped 
and ablated. However, as in the WD-WD case, a^ is expected 
to be in the range ^0.1- 10 AU, and this may preclude a large 
effect on the tertiary, or the SGRB lightcurve. 

5.6. Gravity Waves 

' The fraction of the mass stripped from the ter- 
tiary is ~ lO-3(MsN/M0)(M0/m2)[(R/R0)/(a2/O.lAU)]2 
(vsN/10''km/s)(600km/s/vcsc), where Msn is the mass of the super- 
nova ejecta, i'sn is the ejecta velocity, and I'csc is the escape velocity from 
the tertiary. The mass shock heated and ablated from the tertiary is ~ 6 times 
higher (see Tab. 1 of Wheeler et al. 1975). 



The Diffuse Background: If most close WD-WD binaries 
were born in triple systems, predictions for the gravitational 
wave background may be modified from fiducial predictions 
(e.g.. Farmer & Phinney 2003). For highly eccentric orbits 
the total power in GWs, and the peak frequency of GWs, in- 
creases significantly. The peak GW frequency is well approx- 
imated by (Wen 2003) 

As an example. Figure |6] shows the early (left panel) and late 
(right panel) time evolution of the GW frequency for the triple 
system shown in Figure [T] Although the inner binary, has 
/gw — 0.008 mHz in the absence of a tertiary, with the addi- 
tion of the tertiary few increases strongly, and periodically, 
on a timescale Jk (eq.|4|i. 

As shown in Appendix lAl the maximum eccentricity at- 
tained by the inner binary can be estimated semi-analytically. 
It is these estimates that are used in the calculations of the 
merger time shown in Figures[3]and|5] Substituting this maxi- 
mum eccentricity into equation (|9]l yields /q^, the maximum 
GW frequency produced by the binary during inspiral. The 
horizontal solid line in Figure |6] shows this estimate for /"^^. 

Figure [T] shows a summary of results for computed 
from the results for two panels of Figure[3]for WD-WD bina- 
ries (a2,()/fli,o= 10, 60). For the models that merge in less than 
tu, regions of different peak /^^'^ are shaded in the ai.o-cos/o 
plane. All frequencies shaded are in the range detectable by 
LISA. Note that the WD-WD binaries that attain the highest 
frequencies (interior white) have the highest peak eccentrici- 
ties, the smallest periapses, and will interact tidally, thus mod- 
ifying the results presented here 03.2I ). This can be seen since 
/gw" approaches the inverse dynamical time for an individual 
WD for some regions of parameter space. 

Since the extragalactic GW background is dominated by 
close WD-WD binaries (Farmer & Phinney 2003), which are 
generally assumed to be circular as a result of the preceding 
common envelope evolution, if the fraction of close binaries 
that are actually in triple systems is fairly large, as implied 
by the statistics on solar-type binaries (Tokovinin et al. 2006; 
FT07), then the expected GW background will be modified. 
Recent studies of the LISA GW foreground focus exclusively 
on circular WD-WD binaries (Ruiter et al. 2010), and these 
results too would need to be revisited if the triple fraction of 
close WD-WD binaries is high. The only cases where eccen- 
tric WD-WD binaries have been considered are in globular 
clusters where high eccentricity can be imparted to binaries 
during binary-single and binary-binary interactions (Benac- 
quista 2001; Ivanova et al. 2006; Wiflems et al. 2007). 

Sources & Foreground: Whether NS-NS, WD-WD, or oth- 
erwise, the most eccentric systems will be short-lived, and 
this decreases the probability that they can be seen. The 
competing effect, of course, is that their overall GW lumi- 
nosity is larger The right panel of Figure |6] indicates that a 
fli = 0.05 AU WD-WD binary in a triple system with a 1 Mq 
tertiary at fl2/ai =20 would go through periodic decade-long 
enhancements in ej, that are potentially detectable by LISA 
for nearby systems. During the mission itself, many such bi- 
naries might be seen, both Galactic WD-WD systems with 
larger oi and NS-NS systems with higher /q^- The effect 
of the Kozai mechanism is to make otherwise unobservable 
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Galactic WD-WD binaries observable during pericenter pas- 
sage, since it is at these times at high e\ that is maxi- 
mized. This may cause individual sources to be observable 
at pericenter passage, and it may cause an added degree of 
complexity to the GW foreground from all local sources. See 
Gould (2011). 
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APPENDIX 

APPROXIMATION TO THE MERGER TIME 

I use the methods discussed in MH02 and Wen (2003) to make simple, but accurate, estimates of fmeige to supplement the direct 
calculation of the time-dependence of the orbital elements, as described in ^ Since this scheme provides an efficient way to 
estimate the merger timescale for a very wide range of system parameters, as in Figures [3] & ID I repeat the steps here. 

The goal is to estimate fmerge from the initial conditions of the triple system. Following Wen (2003), I use the fact that 
(neglecting gravitational radiation), the quadrupole-level Hamiltonian is conserved throughout the evolution. It can be written in 
terms of e = l-e^, i, and g\ as (MH02) 

, 9 T 4 (mn + m\) f b2\^ f 2G(mc\ + m\)\ 
H^(e,gi) = -2e+ecos2/ + 5(l-e)sin2gi(cos2/-l) + — -[-] o ]^ (^1) 

where b2 = a2(l—e^^l^, and the last term is the Ist-order post-Newtonian correction that accounts for GR precession. 

Starting with initial values of the eccentricity of the inner binary and its argument of periastron, ei o (eo) ™d g^ Q, as well as 
the initial inclination /q, the system evolves to a maximum eccentricity ei.max and thus minimum e^in^ at a critical ,gi,crit and i^^^. 
Taking de/dt = at the moment e = e^i„ provides a relationship between emin, ,gi.crit, and iait (Wen 2003): 

• ,T s 8 /G(mo + mi)y''VGmomi\ /fl2\ V425-121eminV 2x3/2 .aox 
sm(2gi,erit) = TXT —2 —T— - —r. r- (1-^2) ^ • (A2) 



Using the fact that (Wen 2003) 



225 \ flic^ / \a]chn2 J \ai J \ sin^ lane 



eo \^^^ a 1 

cos icrit = cos /o ( I + — r- , (eo " Cmin) , (A3) 



En 



2^0 VaiCn 



where /3o = ('W2/'Wo'Wi)[('Wo + 'Wi) fl2/(('Wo + 'Wi + 'W2)(l-e2))] ' , allows one to write sin(2gi_crit) (equation IA2l l in terms of just 
emin and the initial parameters of the system. In calculating gicrit and cos/crit, I assume that oi, 02, and 62 are unchanged from their 
initial values. Solving the implicit equation AW = W(€o,gi^o)-W(€mm,g\cnt) = for emin, then allows for an accurate estimate of 
the merger time (MH02): 

-I /2 

fmerge = fGw(ai,0,emin)eniin , (A4) 

where here few is computed from the formalism of Peters (1964). In practice, I solve the equation AW = using Newton-Raphson 
iteration with an initial guess for e^in using the approximate expressions from Wen (2003) and MH02. 

A comparison between the numerically calculated value of Wrge and this approximate scheme is shown in Figure lATI Typically, 
fraei-ge as Calculated from equation (IA4l i is underestimated by a factor of ^ 1.5-2. There are systematic deviations as a function of 
both / an d 02/ a\- These deviations are not due to the fact that the numerical solution is to the octopole-level equations, whereas 
equation (IA4I I is only at quadrupole order, since for mo = nii the equations reduce from octopole to quadrupole and this hypothesis 
can be explicitly tested. The right panel shows that for retrograde orbits and small 02/01 < 20 the estimate of equation ( IA4l i can 
overpredict fmerge by a factor of ~ 10. Note, however, that it is precisely in these cases where tidal effects (neglected here; © are 
expected to be most imp ortan t to the evolution. Future work should address these systematic differences between the fmerge and 
the estimate of equation (IA4b . For the purposes of a quick and fairly accurate estimate of fmerge (as in Figs.[3]and|4]i, Figure lATI 
shows that equation ( IA4b is sufficient. 
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Fig. A1. — Left panel: Comparison between /merge calculated by solving the full system of equations and the estimate of (merge given in equation jA4t as a 
function of a2.o/'^i.o. "l.o = 03 (red), 0.06 (blue), and 0.09 (green), for = 89° (filled circles), 85° (open circles), and 80° (crosses). The estimated merger 
timescale typically undeipredicts ?merge by a factor of ~ 1.5 — 2.5, with systematic offset that depends on ; and ot/"!- Right panel: Same comparison, but as a 
function of cos(i), for a\ = 0.01 AU, and 02/01 = 10 (red), 20 (orange), 30 (yellow), 40 (green), 50 (turquoise), 60 (blue), and 70 (violet). Compare with right 
panel of Fig.|2] 



